Zero temperature phase diagram of the square-shoulder system 
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Particles that interact via a square-shoulder potential, consisting of an impenetrable hard core 
with an adjacent, repulsive, step- like corona, are able to self-organize in a surprisingly rich variety of 
rather unconventional ordered structures. Using optimization strategies that are based on ideas of 
genetic algorithms we encounter, as we systematically increase the pressure, the following archetypes 
of aggregates: low-symmetry cluster and columnar phases, followed by lamellar particle arrange- 
ments, until at high pressure values compact, high-symmetry lattices emerge. These structures are 
characterized in the NPT ensemble as configurations of minimum Gibbs free energy. Based on sim- 
ple considerations, i.e., basically minimizing the number of overlapping coronae while maximizing 
at the same time the density, the sequence of emerging structures can easily be understood. 
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I. INTRODUCTION 



For more than 25 years considerable effort has been 
dedicated to study the thermodynamic, structural, and 
dynamical properties of hard core particles with an adja- 
cent soft repulsive shoulder, i.e., so-called core softened 
potentials. This class of potentials was probably first 
considered by Hemmer and Stell in 1970 [1] in a model 
where the soft repulsion was characterized by a linear 
ramp with an additional, weak attractive tail. The sys- 
tem was introduced in an effort to study the possibil- 
ity of the occurrence of more than one critical point in 
the phase diagram of a simple model system. In numer- 
ous, subsequent investigations evidence was provided for 
a surprisingly rich variety of rather unusual properties of 
this class of systems: these features range from isostruc- 
tural solid-solid transitions, where possibly several solid 
structures are involved (see, e.g., [E 0, M IE ) i_P ver 
a very complex phase diagram of the solid phases [E Q , 
to different sorts of anomalous behavior, encountered in 
the static and/or in the dynamic properties (see, e.g., 

nam in mis)- 

Within this class of core softened potentials the square- 
shoulder interaction, consisting of an impenetrable hard 
core with an adjacent, repulsive shoulder (or corona), 
is undoubtedly the simplest representative. Despite its 
simple, radially symmetric functional form this system is 
nevertheless able to offer a large variety of unexpected 
features, which are mostly related to its structural prop- 
erties. This propensity and ability for unconventional 
self-assembly scenarios was already discovered in the re- 
markable study by Jagla [9] on a particular family of 
iwo-dimensional core softened systems, where the square- 
shoulder system was included as a special case: in this 
contribution evidence was given that the particles are 
able to self-organize in a surprisingly broad variety of 
highly complex ordered structures. In subsequent work 
on the two-dimensional case these particle arrangements 



that include, amongst others, cage- or lane-formation as 
well as micellar or inverse micellar configurations were 
confirmed or newly discovered both in computer simu- 
lations [IE, HE] and in theoretical investigations [ijj • A 
more systematic study of the ordered particle configura- 
tions of the two-dimensional square-shoulder system was 
presented in (l8j . 

The aim of the present contribution is to investigate 
in a systematic and thorough way the ordered particle 
arrangements of the square-shoulder system in three di- 
mensions. To this end we study the system at T = 
and in the NPT ensemble; thus, we search for configu- 
rations that minimize the Gibbs free energy, which we 
will term - to be consistent with previous contributions 
- as minimum energy configurations (MECs). To pro- 
vide a deeper insight into the self-assembly strategies of 
the system we have considered a small, an intermedi- 
ate, and a large shoulder range. While preliminary re- 
sults have already been presented in [IE], we identify in 
the present, more systematic investigation an overwhelm- 
ingly rich variety of MECs. Analyzing these data, we 
give evidence that these MECs can be grouped together 
in four structural archetypes, that emerge in dependence 
of the value of the pressure, P, that is exerted on the 
system: cluster structures are preferentially formed at 
low P- values, while columnar and lamellar structures are 
predominantly identified at intermediate pressure values; 
finally, compact particle configurations emerge at high 
pressure. While this general rule might still be less obvi- 
ous at small shoulder range, it is nearly perfectly obeyed 
for an intermediate shoulder width and definitely holds 
for the case of a broad corona. With its simple func- 
tional form the inter-particle interaction offers not only 
many computational advantages (20| . It allows to un- 
derstand via simple geometrical considerations the sys- 
tem's self-assembly strategy: it is in particular the range 
of the shoulder that turns out to be responsible in a 
highly sophisticated way for the formation of the complex 
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structures. In addition, the flat energetic plateau of the 
shoulder with its finite range represents a very sensitive 
antenna to distinguish between energetically competing 
structures. Our observations provide a deeper insight 
into the system's strategy to form ordered equilibrium 
particle configurations, a knowledge that might be useful 
to understand self-assembly processes in other systems 
with more complex interactions. 

Objections against the simple functional form of the 
potential are refuted by the argument that it is able to 
capture the essential features of colloidal particles with 
core-corona architecture as they are, for instance, treated 
in [2l[ and references therein. And indeed, several of the 
MECs that we could identify had already been encoun- 
tered in previous theoretical, experimental, and com- 
puter simulation investigations: micellar and inverse mi- 
cellar structures Qjl H2] , spirals [23[ , chains and layers 
@, S, HO, El HI H, and cluster phases [H, S3, S3 , to 
name a few examples. 

Although the identification of MECs represents a 'sim- 
ple' optimization problem of the Gibbs free energy, its 
solution has turned out to be highly non-trivial. In this 
contribution we present a systematic sequence of ordered 
MECs for the three selected values of shoulder width. 
This achievement is mainly due to our search strategy, 
which is based on ideas of genetic algorithms (GA). In- 
troduced already several decades ago in a completely dif- 
ferent context [29( | , these approaches have meanwhile be- 
come a highly appreciated optimization tool to identify 
ordered particle arrangements both in hard I3Q . |3U l32fl 
as well as in soft matter systems [H, Ojl |H, Hj, |3f| l36j . 
The high reliability, flexibility, and efficiency of GA-based 
optimization strategies, in combination with a particular 
search strategy that is intimately related to the simple 
functional form of the square-shoulder potential, make 
us believe that the sequences of MECs, that will be pre- 
sented and discussed in the following, are complete. 

The paper is organized as follows. In section |TT] we 
briefly present the square-shoulder system. The subse- 
quent section deals with the theoretical tools of this con- 
tribution: the GA-based search strategy as well as the 
theoretical considerations to identify close-packed parti- 
cle arrangements of the system as a function of the shoul- 
der width. The results of our investigations are summa- 
rized in section IIVI we start with the close-packed par- 
ticle arrangements (as they play a key-role in the search 
strategy) and present and discuss in the following the 
MECs that we have identified for the three different cases 
of shoulder width. The conclusions of the contribution 
are summarized in sectionfVl which also contains the dis- 
cussion of possible future work. 



II. MODEL 

We consider a system of particles that interact via the 
square-shoulder potential, which we parameterize as fol- 



lows: 

{00 r < a 

e a < r < Act . (1) 

Act < r 

a is the diameter of the impenetrable core and Act is 
the width of the adjacent, repulsive shoulder (or corona) 
of height e, e > 0. Further, we introduce the number- 
density, p — N/V (N being the number of particles and 
V being the volume of the system) , and the dimension- 
less number-density, p* = pa 3 . Thermodynamic quanti- 
ties will be used in the following reduced units: pressure, 
P* = Pa 3 /e, internal energy, e* — E/(Ne), and Gibbs 
free energy, g* — G/(Ne). Since we perform our investi- 
gations at T = 0, G = E+PV and hence g* = e*+P*/p*. 

The simple functional form of the square-shoulder po- 
tential with its constant potential barrier and its fi- 
nite r ang e makes the system the 'quintessential test sys- 
tem' [201 ] for the purpose of the present contribution. 
It also simplifies considerably thermodynamic considera- 
tions. For a given periodic particle arrangement, which 
we characterize by the number of overlapping coronae, 
e* is a rational number: it is given as the ratio of the 
number of overlaps per particle in the unit cell, divided 
by the number of these particles, which we denote as rib- 
For this particle arrangement g* = e* + P* / p* , is there- 
fore a linear function of the pressure, P*, and is conse- 
quently represented in the (g*, P*)-plane by a straight 
line: its slope is given by l/p*, while its intercept is the 
energy of the configuration, e*. The limiting particle ar- 
rangement at low pressure is easily identified as a close- 
packed arrangement of spheres with diameter Act: thus 
e* = and the slope of g* in the (g*, P*)-diagram is 
given by l/V^m = A 3 /\/2. For the high pressure limit 
the situation is more delicate: while the slope of g* as 
a function of P* is easily identified for obvious reasons 
as l/p* nax = l/v^, the value of e* depends in a sensi- 
tive way on A. In subsection IIII Bl we will give evidence 
that the square-shoulder system shows a rich variety of 
close-packed scenarios as A varies. 

With the above considerations in mind, we can antic- 
ipate that the $*-values of all MECs will be located on 
a sequence of intersecting straight lines in the (g*,P*)- 
plane, each of them being characterized by a slope of 
l/p*, with 1/ '/Omi n > > 1/Pmax, an d by an intercept 
e*. This fact will simplify considerably our search for 
MECs (see discussion in subsection II V Bp . 

III. THEORY 

A. Genetic algorithms 

The MECs of our system have been identified with a 
search strategy that is based on ideas of genetic algo- 
rithms (GAs). GAs are very general optimization tools 
that model natural evolution processes, such as recom- 
bination, mutation, or survival of the fittest (29| . Their 
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successful applications in a wide range of fields demon- 
strate their flexibility and reliability. The basic ideas of 
GAs can be summarized as follows: the central quantity 
of this concept is an individual, Z, which represents a pos- 
sible solution to the problem. The quality of a solution, 
i.e., of an individual T, is measured via a so-called fitness 
function, f(T). Individuals with a higher fitness value 
are assumed to be of higher quality. In our search for 
ordered particle configurations that minimize the Gibbs 
free energy, an individual corresponds to a lattice while 
the fitness function is related to G and will be speci- 
fied below. Starting from a large number of individu- 
als, which represent the initial generation, individuals of 
a subsequent generation are created with recombination 
and mutation processes, both of them having a highly 
stochastic character. Individuals with a higher fitness 
value are preferred in the reproduction process. In addi- 
tion, mutation operations are performed on the individ- 
uals with some probability p m . By iterating this process 
we create a reasonably large number of generations. The 
final result of the GA-based search strategy is the indi- 
vidual with the overall highest fitness value. 

For our particular problem an individual T is identi- 
fied by a (possibly non-simple) periodic crystal structure. 
Due to the highly stochastic character of the reproduc- 
tion and of the mutation processes, a straightforward im- 
plementation of the algorithm is prone to propose a large 
number of ordered configurations where the hard cores of 
the particles overlap and which therefore correspond to 
unphysical particle arrangements. This, in turn, causes a 
drastic reduction of the efficiency of the algorithms. To 
overcome this problem we have developed a particular 
parameterization of an arbitrary simple lattice via three 
lattice vectors {ai,a2,a3} [371 ] : here, a\ — |ai| repre- 
sents the shortest possible distance between two lattice 
sites in the entire lattice and a 2 — |&2 1 is the second 
smallest distance in the lattice (i.e., a\ < 02) with ai 
and a2 being linearly independent; finally, a similar rela- 
tion holds between a3 on one side and ai and a2 on the 
other side. Thus if a\ > a, it is guaranteed that the hard 
cores of the particles will not overlap and consequently 
the GA will create only simple lattices where overlap of 
the cores is avoided a priori. For non simple lattices, the 
distances between all particles within the unit cell and 
including also the particles of the 26 neighboring cells 
have to be determined. If the smallest of these distances, 
ioi is smaller than 01, then the lattice is scaled with a 
factor ai/Zo. The rather complex formalism is most con- 
veniently implemented in the NPT ensemble. Thus, a 
state is characterized by a value for the pressure, P, while 
the equilibrium density, p, is a result of the optimization 
procedure. For details we refer to [37| . 

For the implementation of the individuals we use the 
encoding strategy presented in [35[ . Distances and angles 
are encoded in binary strings with a length ranging from 
four to six. Since in the NPT ensemble the Gibbs free 
energy has to be minimized, we use the following fitness 
function 



f(l) = exp{-[G(Z) - G(2b)]/G(2b)}, (2) 

where To corresponds to some reference structure. A pool 
of 700 individuals is evolved through reproduction and 
mutation processes over 500 generations: creation of indi- 
viduals of the new generation from individuals of the pre- 
ceding generation is carried out via one-point or random 
cross-over operations while the mutation process, which 
re-introduces lost genetic materials and avoids inbreed- 
ing, is realized with a mutation probability of p m = 5%. 
For each state point, 1000 of such independent runs have 
been carried out. Finally, the individual with the over- 
all lowest G-value, X m j n , is considered to be the solution 
of the GA. To account for the limited accuracy caused 
by the encoding procedure, the parameters of 2"mjn were 
refined via a final Powell optimization algorithm [38j |. 

B. Close- packed structures 

The limiting case at high pressure is always a lattice 
where the hard cores of the particles arrange in a close- 
packed structure. Thus the slope of the line that ex- 
presses the linear dependence of g* on P* in the (<?*, P*)- 
plane is obviously given by l/p^ax — l/v2- However, 
the intercept of this line, i.e., the energy of this arrange- 
ment, e*, requires more careful considerations. Below 
we will give evidence, that the square-shoulder system 
is able to self-organize not only in the well-known close- 
packed scenarios, i.e., in fee or hep lattices, but also in 
more complex structures [3!|- We emphasize that the 
square-shoulder system serves - due to its flat plateau 
and due to the finite range of the corona - as an antenna 
that is able to identify in a very sensitive way between 
competing particle arrangements. 

In an effort to find for a given value of A the ener- 
getically most favorable close-packed arrangement of the 
particles we proceed as follows. We consider the lattice 
as being built up by periodically repeated stacking se- 
quences of ni hexagonally close-packed layers, introduc- 
ing for convenience the conventional labels, A, £?, and C 
[40] . A stacking sequence of ni layers can therefore be 
described by a string of ni of these symbols. The trivial 
close-packed arrangements, fee and hep, arc thus char- 
acterized by the sequences ABC (with ni = 3) and AB 
(with ni = 2). For a given value of ni we consider all 
possible stacking sequences of length m ; without loss of 
generality we start all sequences with the label A. Some 
of the proposed sequences have to be ruled out: this is 
the case when two neighboring layers carry the same in- 
dex. Some of them can be ruled out: this is, for instance, 
the case when symmetry considerations reveal that two 
different stacking sequences lead to the same lattice. 

Pursuing this strategy we find for the smallest up- 
values the following situation: for ni — 2 we have only 
the hep structure (AB) and for ni — 3 we recover the 
fee lattice (ABC). Also for four- and five-layer stackings 
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only one representative remains: ABAC and ABABC 
can be identified, respectively. At n./ = 6, we encounter 
for the first time two non-equivalent stacking sequences, 
namely ABABAC and ABACBC. With increasing n/ 
the number of possible stacking sequences increases dras- 
tically. For instance, for A = 4.5, where we have con- 
sidered stackings with up to 12 layers, we were able to 
identify 136 different sequences. A comprehensive table 
of possible stacking sequences for a given value of ni is 
presented in (4lT |. 

Finally, for a given value of A we include a sufficiently 
large number of layers and evaluate and compare the en- 
ergies e* of all candidate stackings. The finite range of 
the shoulder and its flat energy plateau help to reduce 
the numerical effort considerably. On most occasions we 
encounter degeneracy, i.e., two (or even more) different 
stacking sequences are characterized by exactly the same 
value of e*. In such cases, we consider the shortest among 
these stacking sequences to be the energetically most fa- 
vorable configuration with the only exception that we 
favor fee to hep [42j]. 
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FIG. 1: Fig. [T] Energy per particle, e* , for the simplest, en- 
ergetically most favorable close-packed particle arrangements 
for the square-shoulder system as a function of A (full curve). 
Symbols specify the stacking sequences as labeled; see also 
text. 



IV. RESULTS 

A. Close-packed structures 

With the above considerations in mind we can now 
identify the equilibrium close-packed structures for the 
square-shoulder system as they occur at high pressure 
values. These particle arrangements are summarized in 
Figure [T] starting from A = 1 (corresponding to hard 
spheres) and extending to a shoulder width of A = 4.5. 
The figure contains the energy e* of the respective struc- 
tures and symbols characterize their stacking sequences. 

As long as 1 < A < \/2, only nearest neighbor inter- 
actions have to be considered. All stackings are charac- 
terized by the same number of overlapping coronae; the 
energy of a tagged particles amounts to half the num- 
ber of nearest neighbors, i.e., e* = 6. Although for 
\J~2 < A < 2-^/2/3, the second nearest neighbors start to 
play a role, we still obtain for all possibl e sta ckings the 
same value for e*, namely e* — 9. For 2-^/2/3 < A < \/3, 
the different stackings are characterized by different en- 
ergy values; among these, the fee structure is one of 
the stackings with the lowest e*-value, namely e* = 9, 
while all other stackings are energetically equal or less 
favorable. Therefore in this A-interval, fee remains the 
simplest, energetically most favorable structure. For 
A = y/3, a hep lattice with e* = 19 becomes the sim- 
plest structure with the lowest e-value. Further stacking 
sequences can be extracted from Figure [1] for A- values up 
to 4.5. This figure shows that for a few A-intervals also 
stacking sequences other than fee or hep are obtained 
as the energetically most favorable close-packed particle 
arrangements at high pressure. 



B. Configurations that minimize the Gibbs free 
energy 

Once we have determined the limiting high pressure 
MECs, we can proceed to the identification of the whole 
sequence of MECs as a function of the pressure. This 
is done in the following three subsections where we have 
considered square-shoulder systems with a short (A = 
1.5), an intermediate (A = 4.5), and a large (A = 10) 
shoulder width. Abbreviations of the underlying lattices 
that are used in the text and in the Figures are summa- 
rized in Table |TJ 

For A = 1.5 we shall give a detailed geometrical inter- 
pretation of these particle arrangements; this will pro- 
vide clear evidence about the system's strategy to ar- 
range the particles at given pressure in such a way as 
to minimize the number of overlapping shoulders and to 
maximize at the same time the density. Although we 
will not be able to pursue these geometrical considera- 
tions in full detail for the other A-values, we will be able 
to identify an emerging sequence of structural archetypes 
as we increase the pressure: while at low pressure parti- 
cles tend to arrange in clusters, which then populate the 
positions of regular lattices, we encounter with increas- 
ing pressure columnar, lamellar, and, finally, compact 
structures. With a few exceptions this rule is obeyed for 
A = 4.5, while it is strictly followed for A — 10. For a 
more detailed presentation we refer to [4l| . 

Before we present and discuss the sequences of MECs 
in detail we briefly outline how we can take benefit from 
the fact that for this particular system g* is a linear func- 
tion of P*. Our search algorithm is sketched in Figure 
[2] In a first step, we determine the intersection point of 
the two straight lines in the (g*, P*)-diagram which rep- 
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FIG. 2: Schematic representation of our search strategy to 
identify MECs in the (g*, P*)-plane. The dotted lines rep- 
resent g* as a function of P* for the limiting low and high 
pressure configurations: The vertical arrows represent GA- 
runs that identify, starting from an initial guess (dot), an 
energetically more favorable MEC (square). For details see 
text. 



fident that the sequences of MECs that we shall present 
in the following, are complete. 



1. Short shoulder width (X = 1.5) 

The phase diagram (i.e., g* and e* as functions of P*) 
for the square-shoulder system with A = 1.5 is depicted 
in Figure O the corresponding ordered equilibrium struc- 
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resent the high and the low pressure limiting cases; let 
the corresponding pressure value be Pr*i. At this state 
point we perform a sequence of GA searches. In each of 
these runs we consider a different number of basis atoms, 
where the maximum number of basis particles depends 
on the value of A. This optimization step leads to a new 
particle configuration which is characterized by a Gibbs 
free energy that is lower than the one of the intersec- 
tion point, by a density p*^, and by an energy e*^. Thus 
this particle arrangement is at given pressure P^i the en- 
ergetically most favorable one. ef^ and p*^ define a new 
line in the (g*, P*)-plane; we determine the two inter- 
section points of this line with the two lines representing 
the limiting configurations leading to the pressure val- 
ues Pj2 a i and P&b] . At these two state points we launch 
new GA searches. This procedure is repeated until at 
none of the intersection points of an iteration step an 
energetically more favorable particle arrangement can be 
identified. On one side this procedure avoids a rather 
time-consuming scan of the pressure range on a finite 
grid and thus brings along a considerable reduction of 
the number of GA steps and, consequently, of the com- 
putational effort; on the other side this strategy avoids 
the risk of simply 'forgetting' MECs. Both features be- 
come more and more important with increasing shoulder 
width, since the distribution of MECs over the whole 
pressure range is highly nonlinear, as can be seen in Fig- 
ures 131 El and [7] 

This systematic search strategy, in combination with 
the reliability of GA-based optimizations, make us con- 



FIG. 3: g* and e* as functions of P* for the ordered equilib- 
rium structures identified for the square-shoulder system with 
A = 1.5, as labeled. The dotted lines indicate the low and high 
pressure limiting configurations (see text). The identified lat- 
tices are indicated by standard abbreviations (see table HJ , 
including, if required, the number of basis particles; see also 
figure [4] 

tures are compiled in Figure [H except for the trivial 
low and high pressure structures, where particles arrange 
in any close-packed crystal structure, which we take to 
be fee (see discussion above). Further numerical details 
about the seven identified MECs are compiled in Table 
HT1 Although we have considered in our search strategy 
lattices with up to eight basis particles, only lattices with 
at most two basis particles were identified. The limiting 
low pressure configuration is characterized by e* = and 
p* = \/2/X 3 ~ 0.629. Further, since A = 1.5 is slightly 
larger than y/2, we are above the threshold value (cf. dis- 
cussion in subsection IIV Ap where the hard cores of the 
particles form a close-packed structure and only the coro- 
nae of nearest neighboring particles overlap; thus e* — 9 
and p* — 

As we start our search in the low-pressure regime, 
the first non-trivial structure we encounter is a body 
centered orthorhombic (bco) structure [Figure BJl)]. A 
more detailed consideration identifies this particle ar- 
rangement as a columnar structure: particles form lanes 
along which the hard cores are in direct contact. While, 
of course, these lanes lead to an intra-columnar shoulder 
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(1) bco, e* = 1 (2) fco, e* = 2 




(3a) hex, e* = 4 (3b) hex (rect) 




(4a) sfco, e* = 5 (4b) sfco (rect) 




(5a) bco, e* = 7 (5b) bco (rect) 



FIG. 4: Visualization of the non-trivial ordered equilibrium 
structures for the square-shoulder system with A = 1.5. Struc- 
tures are characterized by standard abbreviations (see table 
HI and their respective e*-value. Color code: green - parti- 
cles at the corner positions of the conventional unit cell; red - 
particles at body- or face-centered positions; blue - additional 
basis particles. 



overlap along the lanes, any other inter-columnar over- 
lap is avoided; consequently, e* = 1. Simple geometric 
considerations reveal that the edge lengths of the con- 
ventional bco unit cell have the following values: ct, Act, 
and V3A 2 - 1 ct. 

The avoidance of inter-columnar shoulder overlap has 
to be sacrificed as the pressure is further increased, lead- 
ing to a rather compact structure: we identify a face cen- 
tered orthorhombic (fco) lattice with an additional basis 
particle (color code: blue) [cf. Figure [1(2)]. This par- 
ticle is in direct hard core contact with its four nearest 
neighbors: three of them are located at the faces (color 
code: red) and one sits at the corner (color code: green) 
of the conventional unit cell. Furthermore, this parti- 
cle is separated by a distance Act from its second nearest 
neighbor (color code: green) which occupies another cor- 
ner of the conventional orthorhombic unit cell. Finally, 
particles at the smallest faces of the cell (color code: red) 
are positioned in such a way that their shoulders touch 
the shoulders of the particles located at the corners of 



the corresponding face of the conventional orthorhom- 
bic unit cell (color code: green). These considerations 
fix the edge lengths of the unit cell to be \/2A 2 — 2 ct, 
V2A 2 + 2 ct, and 2%/4 - A 2 ct. 

As we further increase the pressure the particles ar- 
range in a MEC that can be identified as a lamellar struc- 
ture, cf. Figure|3K3a). It can be described as a stacking of 
hexagonally close-packed layers which are placed exactly 
on top of each other. The nearest neighbor distance is 
obviously ct. The second nearest neighbors are separated 
by the inter-layer distance which is fixed by the require- 
ment that the corona of a tagged particle touches the 
shoulders of its twelve third nearest neighbors, located 
in the adjacent layers. Thus, this distance amounts to 
VA 2 — 1 ct ~ 1.12ct. It should be pointed out that par- 
ticles in nearest and second nearest neighbor distance 
form a rectangular particle arrangement [emphasized in 
Figure [4|3b)] that will also be encountered in the subse- 
quent MECs: if we consider within a layer two particles 
in close contact, then they form with the corresponding 
particles of one of the adjacent layers a rectangle with 
edge lengths ct and VA 2 — 1 ct. 

For even higher pressure values only compact struc- 
tures are identified. The next MEC can be described as 
a single face centered orthorhombic (sfco) lattice with 
two basis particles, visualized in Figure E(4a). The 
orthorhombic unit cell is built up by two side faces 
that have exactly the aforementioned rectangular shape 
[formed by particles in green, emphasized in Figure 
0|4b)], while the two larger side faces are each decorated 
in their center by an additional particle (color code: red). 
Finally, the additional basis particles (color code: blue) 
are located in such a way, that they are in direct contact 
both with the four particles forming the side faces as well 
as with the two red particles located in the other side 
faces. Simple geometric considerations lead to the edge 
lengths of the orthorhombic cell, namely: ct, VA 2 — 1 ct, 
and (V3 + V4- A 2 ) ct. 

The last non-trivial compact structure is a body cen- 
tered orthorhombic lattice [emphasized in Figure [U(5a)]. 
Again, we can easily identify the side faces of the conven- 
tional unit cell as the above mentioned rectangular struc- 
ture [see Figure |U(5b)]. In addition, the central particle 
is in direct hard core contact with the particles forming 
the unit cell. Based on these geometrical considerations 
the edge lengths of the unit cell can easily be identified 
to be ct, VA 2 — 1 ct, and y/4 — A 2 ct. 



2. Intermediate shoulder width (X = 4.5) 

A much larger diversity in the ordered equilibrium 
structures could be identified for an intermediate shoul- 
der width of A = 4.5. The limiting low pressure MEC is 
of course again an fee structure with a nearest neighbor 
distance Act and, hence, e* = 0. On the other hand, the 
high-density limiting particle configuration is an hep lat- 
tice with e* = 263 (see Figured] and discussion in subsec- 
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tion llV A| . With the help of the GA we have obtained in 
total 33 different MECs over the entire pressure regime. 
In our investigations unit cells with up to 10 basis parti- 
cles have been considered, in the end only configurations 
with up to 8 basis particles were part of the MECs. 

A first look at these MECs gives evidence that within 
this sequence of MECs we can easily identify the afore- 
mentioned four structural archetypes: at low pressure, 
the system prefers to form cluster structures; with in- 
creasing pressure, columnar structures are formed, which 
then transform into lamellar particle arrangements; fi- 
nally, at high pressure, we observe compact structures. 
This rule, which is disobeyed only twice for A = 4.5, 
can nicely be understood via a detailed analysis of the 
particle arrangements, reflecting the system's strategy to 
reduce at a given pressure the number of overlapping 
coronae as much as possible (i.e., minimizing e*) while 
maximizing at the same time the particle density. 

Numerical details about these ordered structures are 
summarized in Table IIII1 The phase diagram for the 
square-shoulder system with A = 4.5 is depicted in Fig- 
ure [5] It also contains information to which class of the 
four archetypes a particular MEC belongs. Finally, the 
horizontal bar at the bottom of the figure indicates those 
MECs where no direct contact between the cores of the 
particles occurs. 
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FIG. 5: g* and e* as functions of P* on a double logarithmic 
scale for the ordered equilibrium structure identified for the 
square-shoulder system with A = 4.5, as labeled; note that 
due to the non-linear scale, the linear dependence between 
g* and P* is no longer visible. The structural archetypes 
to which a given MEC belongs (cluster, columnar, lamellar, 
or compact structure) are specified by a symbol (as labeled). 
The black horizontal bar at the bottom of the figure indicates 
those ordered particle arrangements where no direct contact 
between the cores of the particles is observed. 

At low pressure- values (i.e., up to P* ~ 1.70) particles 
prefer to arrange in ordered clusters of up to eight par- 
ticles, which, in turn, populate the positions of crystal 



lattices. A closer analysis of these structures reveals a 
strong interplay between the shape of the clusters and 
of the symmetry of the unit cell: the more aspherical 
the clusters are the lower is the symmetry of the lattice. 
This tendency reflects the system's efforts to avoid to a 
highest possible degree a shoulder overlap of neighboring 
clusters. A nice visualization of this strategy can, for in- 
stance, be observed in the structure depicted in Figure 
Et5): the rather elongated four-particle clusters are lo- 
cated on a low-symmetry triclinic lattice; on the other 
hand, the nearly spherically shaped eight-particle clus- 
ters of the structure depicted in Figure populate 
the lattice positions of a bco structure, that has a con- 
siderably higher symmetry. A systematic, quantitative 
analysis of all cluster structures reveals that for most of 
these MECs only rarely shoulder overlaps of neighboring 
clusters are observed. A nice example that demonstrates 
the complexity of cluster structures is depicted in Figure 
(Hl(4). In this MEC we can identify two different sorts of 
clusters: tetrahedral clusters occupy the corners of a tri- 
clinic lattice, while the other four-particle cluster species 
populates a central position in the body of the triclinic 
cell. 

As the pressure is further increased there is a dras- 
tic change in the system's strategy to arrange particles, 
namely the formation of columnar structures, where par- 
ticles self-organize in lanes. This leads to a considerable 
energetic penalty, since - due to the short inter-particle 
distance within the columns - an appreciable number 
of overlapping shoulders is induced; at the same time 
a rather high particle density is guaranteed along these 
lanes which, in turn, contributes to e reduction in the 
Gibbs free energy. Simultaneously, the system tries to 
compensate for this high energetic cost within the lanes 
by arranging these columns in such a way as to mini- 
mize the inter-columnar shoulder overlap. This strategy 
leads first to the formation of single-columnar structures 
[as depicted in Figure [6j7)], and later, as the pressure 
is increased, to double-columnar particle arrangements 
[cf. Figure [SJ 9-11)]. We point out that within the lanes 
particles are only in direct contact at sufficiently high 
pressure; in double-columnar structures, particles of ad- 
jacent columns are always in direct contact. The system's 
strategy to avoid shoulder overlap between the lanes can 
nicely be traced in a closer analysis: for the structures 
depicted in Figures [DJ7,9, 10), no corona overlap between 
the single or double strands is observed; only at suffi- 
ciently high pressure - cf. Figure 01 1) - the coronae of 
different neighboring double columns start to overlap. 

For pressures values above P* > 3, the system has 
to search for new ideas of how to minimize the Gibbs 
free energy. Now the change to a new strategy is con- 
siderably smoother than the preceding one: in an effort 
to cope with the increasing pressure the system forms 
lamellar structures - cf. Figures Et 12-21) and Figure 
[SJ24). These MECs emerge from columnar structures as 
the columns approach each other, forming thereby lamel- 
lae; some intermediate stages of this transition can be 
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FIG. 6: Visualization of all 33 ordered equilibrium structures for the square-shoulder system with A = 4.5. Structures are 
characterized by standard abbreviations (see table UJ) and their respective e*-value. Color code: green - particles at the corner 
positions of the conventional unit cell; red - particles at body- or face-centered positions; blue - additional basis particles. The 
shoulders of the yellow particles in panel (17) touch the ones of the other yellow particles, located in the neighboring layers. 
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observed in Figure [6]( 12) and Figure [6^15). Within the 
lamellar structure the system's strategy is obvious. First 
optimize the packing inside a layer, leading to hexagonal 
particle arrangements inside a lamella: while at low pres- 
sure values [Figures [6^13, 14)] particles are more loosely 
packed, they are forced to form a nearly hexagonally 
close-packed structure with a nearest neighbor distance 
of ~ 1.03er at higher pressure [Figures [5)^16, 17)]. Particu- 
lar attention should be dedicated to the latter structure: 
the three neighboring, parallel planes depicted in Figure 
EJ17) are not equally spaced; the two different emerg- 
ing inter-lamellar distances are rather governed by the 
fact that the shoulders of the particles marked in yellow 
located in the three neighboring layers touch. If the pos- 
sibility for optimizing the packing within a single lamella 
has been exhausted, the system starts to form double 
layers [cf. Figures[6jl8-21)] or even triple layers [cf. Fig- 
ure[Sf24)]. A closer analysis of the double-layer structure 
reveals a very complex strategy which we try to explain 
as follows. We consider two pairs of neighboring double- 
layers. On one hand we observe shoulder overlap between 
single layers (belonging to different pairs) that face each 
other: for instance, in Figure [SJ19), the layer formed by 
blue particles in the left-most layer pair and the layer 
formed by red particles in the central layer pair; on the 
other hand the distance between pairs of layers is chosen 
in such a way as to avoid shoulder overlap of single lay- 
ers that do not face each other: for instance, in Figure 
[^K 19) , the layers formed by green particles and the layers 
formed by red particles. These observations turn out to 
be valid for all double layer structures that have been 
identified for this particular shoulder width. 

Finally, we enter for high pressure values the regime of 
compact structures, characterized, in general, by a large 
number of nearest neighbors. In most of these MECs di- 
rect core contact is avoided (cf. horizontal bar in Figure 
[5]), only in the high- pressure regime, where a ct-lattice 
[cf. Figure EJ32)] and, finally, an hep-structure [the lim- 
iting case for A = 4.5, see Figure [61(33)] are the respec- 
tive MECs, the cores are in direct contact. Again, with 
simple geometric considerations, the system's strategy to 
form MECs can be traced back to avoiding unnecessary 
shoulder overlap while maximizing the particle density. 



tigations and in computer simulations for the PSM [45j 
and via purely theoretical considerations, combined with 
computer simulations for a closely related model poten- 
tial [U [H, [H, H3 . As we will show below, such cluster 
phases can also be observed for the square-shoulder sys- 
tem at low pressure values where the hard cores of the 
particles still have a negligible effect on the properties of 
the system. 

Since the MECs are expected to be rather complex we 
have considered up to 29 basis particles in our GA based 
search strategy, up to 22 appeared in the MECs. In total 
we have identified as much as 47 MECs, i.e., a relatively 
large number which makes a detailed discussion and in- 
terpretation of the structures impossible. The thermo- 
dynamic properties of all these MECs are displayed in 
Figure [JJ We point out that the high pressure limiting 
configuration is a hep lattice with e* = 2947. For the case 
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FIG. 7: g* and e* as functions of P* on a double logarithmic 
scale for the ordered equilibrium structures identified for the 
square-shoulder system with A = 10, as labeled; note that 
due to the non-linear scale, the linear dependence between 
g* and P* is no longer visible. The structural archetypes to 
which a given MEC belongs (cluster, columnar, lamellar, or 
compact structure) are specified by a symbol (as labeled). Or- 
dered equilibrium structures marked by crosses are visualized 
in Figure [8] 



3. Large shoulder width (X = 10) 

Finally, we consider the case of a large shoulder width 
for which we have chosen A = 10. Now the hard core re- 
gion is relatively small with respect to the shoulder range. 
Thus at low densities the core plays a minor role and the 
system becomes closely related to the penetrable sphere 
model (PSM) [43j. The PSM belongs to a class of soft 
matter systems where particles are able to solidify in so 
called cluster phases j44j, i.e., where stable clusters of 
particles form, which, in turn, populate the positions of 
periodic lattices. Evidence for this particular phase be- 
havior has been given in density functional based inves- 



A = 10 the rule for the sequence of structural archetypes 
(cluster - columnar - lamellar - compact structures) is 
strictly obeyed (see symbols in Figure [7]). 

As expected cluster structures emerge at low pressure 
values. A few examples of the ten cluster structures that 
have been identified are depicted in Figures [SJ 1-3). The 
clusters can contain as many as 22 particles [e.g., in the 
structure in depicted in Figure [5^3)] and are arranged 
in complex structures. An example for a typical cluster 
is depicted in Figure [8^2); in general the intra-cluster 
arrangement of the particles turns out to be irregular. 

At P* ~ 0.76 the transition to the columnar structures 
occurs. The relatively large shoulder width allows for a 
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FIG. 8: Visualization of a selection of the 47 ordered equilibrium structures for the square-shoulder system with A = 10. 
Structures are characterized by standard abbreviations (see table[I|) and their respective e*-value. Color code: green - particles 
at the corner positions of the conventional unit cell; red - particles at body- or face-centered positions; blue - additional basis 
particles. 



large variety of columnar morphologies, including multi- 
columnar arrangements or complex helical columns - cf. 
Figures [U 4-6). In Figure[8][5) a side view of a single col- 
umn gives evidence of its complex internal structure. Ten 
basis particles were required to parameterize this MEC; 
a single column can be considered to be built up by a 
sequence of aligned clusters. We point out that also in 
experiment helical columns were have been observed for 
a particular class of colloidal particles [23] . Examples for 
multi-columnar arrangements are the triple-columns dis- 
played in Figure H{4), or the MEC shown in Figure [5f6): 
here six parallel single-columns that are nearly in close 
contact are aligned in parallel to build the sixfold-column, 



as can be seen from the right-most column, where the di- 
rection of projection has been chosen to be parallel to the 
columnar axis. 

Most of the MECs identified for the case A = 10 have 
lamellar character: in total we have identified as much as 
28 lamellar MECs. Again, we observe a similar strategy 
as the one identified for A = 4.5: first the particle ar- 
rangement within the single layer structures is optimized; 
then, if this possibility for close-packed arrangements is 
exhausted, multi-layer structures are formed. The large 
shoulder width is responsible both for the large inter- 
layer distance as well as the close contact within groups of 
lamellae: One has the impression that the large range of 
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the shoulder compactifies adjacent layers, bringing them 
in direct contact, while maximizing at the same time the 
distance between these groups of layers [cf. Figures Ej8- 
12)] ; 

Finally, we enter the regime of compact structures. 
Since they resemble very closely those MECs that have 
been identified for A = 4.5, we do not present them here. 



V. CONCLUSIONS 

In this contribution we have thoroughly investigated 
the phase diagram of the square-shoulder system at T = 
0, taking into account a short, an intermediate, and a 
large shoulder width. Measuring the range of the corona 
in terms of Act (where a is the hard core diameter) we 
have assumed the following specific values for the three 
cases: A = 1.5, A = 4.5, and A = 10. Investigating 
the system in the NPT ensemble we have searched for 
ordered particle configurations that minimize the Gibbs 
free energy; this means that the internal energy is mini- 
mized while the particle density is simultaneously maxi- 
mized. These particle arrangements have been identified 
by means of a search strategy that is based on ideas of ge- 
netic algorithms. With this reliable, flexible, and efficient 
optimization tool at hand and taking benefit of the simple 
functional form of the inter-particle potential, which con- 
siderably facilitates both numerical calculations as well as 
geometrical interpretations, we give evidence that the se- 
quences of emerging particle configurations of minimum 
energy are complete. 

A first look on the total of configurations gives clear 
indications that the formation of the particle arrange- 
ments follows well-defined rules as the pressure is in- 
creased: while at low pressure values the system prefers 
to form clusters, which, in turn, populate the positions of 
low-symmetry lattices, we encounter at medium pressure 
values columnar and then lamellar structures. As a rule 
of the thumb we found, the distances between clusters, 
columns and lamellae are always roughly equal to the 
shoulder width. Finally, at high pressure values, rather 
compact particle configurations are identified, that are in 
general characterized by a large number of nearest neigh- 
bors. While this rule for the structural ordering might be 
still less obvious for small values of A, it becomes more 
apparent with increasing range of the corona: at the in- 
termediate A-value of 4.5, it is disobeyed only at two 
occasions, and for A = 10 the ordered particle configura- 
tions fully match this rule. 

The large variety of ordered equilibrium configurations 
that have been identified is overwhelming. It represents 
an impressive example of the capacity and propensity of 
soft matter particles to self-organize in highly non-trivial 
structures. For demonstration we pick out two particular 
examples: first a cluster structure, where the clusters are 
composed of 22 particles, which, in turn, populate the po- 
sitions of a triclinic lattice; second a columnar structure, 
where the particles align in a complex, helical column. 



The particular shape of the repulsive shoulder, i.e., 
its flat energetic plateau in combination with the well- 
defined range of the corona, provides in addition the 
unique possibility to understand the system's strategy 
to form these complex structures. A detailed structural 
and energetical analysis of the emerging configurations 
reveals that the shoulder width plays the dominant role 
in this process; however, also the hard core can have 
considerable influence on the structure formation since it 
represents a lower boundary to inter-particle distances. 
Thus, this study can be viewed as a pedagogical exam- 
ple which provides a deeper insight why particles arrange 
at a given state point in a particular structure, a knowl- 
edge, which might be helpful in the investigations of self- 
assembly processes of systems with more complex inter- 
particle potentials. 

What are the next steps? The most obvious extension 
should be directed towards the investigations of the phase 
diagram at finite temperature, i.e., one should address 
the question which of these configurations will 'survive' 
at T > 0. The theoretical route to this answer is quite 
straightforward: one has to 'simply' merge the proposed 
search strategy with a suitable method to evaluate the 
thermodynamic properties of the system at a finite tem- 
perature. The obvious candidate to evaluate the ther- 
modynamic properties of the sy_stem is of course classi- 
cal density functional theory [481 ] . However, we have to 
raise immediately two serious concerns. First, no reliable 
density functional format for the square-shoulder system 
is available at present; treating, alternatively, the hard 
core within a suitable fundamental measure theory for- 
mat (as, e.g., the one proposed in (49j |) and considering 
the shoulder by a mean-field type perturbative approach, 
risks to provide data that are not sufficiently accurate 
to reliably distinguish between energetically competitive 
structures. Second, in search strategies based on ideas of 
genetic algorithms, the evaluation of the fitness function 
for an individual represents the numerical bottleneck in 
this approach. According to our experience, a combina- 
tion of classical density functional theory with this par- 
ticular search strategy might still be too time-consuming 
for present-day computers. Alternatively to the theoret- 
ical route one might of course perform computer simu- 
lations and we point out that recently investigations in 
this direction have been performed on a related system 
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VI. TABLES 



TABLE I: Standard abbreviations for the 14 Bravais lattices 
used in the text and the captions. The seven crystal systems 
are separated by double lines. 





Oi \j Ul L V IdLlUll 


simple cubic 


SC 


UUUy L.l^llLL.l LAI L^UUlL-< 


bee 


face centered cubic 


fee 


hexagonal 


hex 


trigonal (rhombohcdral) 


trig 


simple tetragonal 


St 


centered tetragonal 


ct 


simple orthorhombic 


SO 


single face centered orthorhombic 


sfco 


body centered orthorhombic 


bco 


face centered orthorhombic 


fco 


simple monoclinic 


sm 


single face centered monoclinic 


sfcm 


triclinic 


trie 



TABLE II: Numerical details of the ordered equilibrium struc- 
tures identified for the square-shoulder system with A = 1.5: 
the underlying lattice is characterized by the according ab- 
breviation (cf. Table [TJ) , m is the number of basis particles 
required to describe the MEC. e* and p* are the energy per 
particle and particle density, respectively. Since the MECs 
can be interpreted on the basis of geometric considerations, 
p* can be given in closed, analytic expressions. 



lattice 


n b 


e* 


P* 


fee 


1 





^ ~ 0.419 


bco 


1 


1 


8 _ ~ 0.556 

3V23 


fco 


2 


2 




hex 


1 


4 




sfco 


2 


5 


16 ~ 1 17 

2 VTE+V35 


bco 


1 


7 




fee 


1 


9 


V2 ~ 1.41 



TABLE III: Numerical details of the ordered equilibrium 
structures identified for the square-shoulder system with A = 
4.5: the underlying lattice is characterized by the according 
abbreviation (cf. Table U}, rib is the number of basis particles 
required to describe the MEC. e* and p* are the energy per 
particle and particle density, respectively. The abbreviation 
in the third row indicates to which of the four archetypes the 
MEC belongs (cluster, columnar, lamellar, compact). Closed, 
algebraic expressions for p* can be derived, but are not pre- 
sented here due to space limitations. For details cf. [4lj . 



lattice 


rib 


shape 


e* 


P* 


fee 


1 


clu 





0.0155 


sm 


2 


clu 


1/2 


0.0235 


sfcm 


3 


clu 


1 


0.0291 


trie 


8 


clu (4) 


3/2 


0.0343 


trie 


4 


clu 


7/4 


0.0369 


trie 


6 


clu 


5/2 


0.0454 


bco 


1 


col 


3 


0.0512 


bco 


8 


clu 


35/8 


0.0638 


trie 


2 


col 


13/2 


0.0799 


Gill 


2 


col 


8 

O 


0923 


LI l\-> 


2 




21/2 


1 04 


bco 


1 
1 


lam / coi 


10 


u. 100 


trig 


1 


lam 


18 


U.104 


trie 


z 


lam 


on 


U.lDo 


trie 




£i 


1 P TYl I fd \ 1 


4Q 10 




trig 


2 


lam 


30 


0.243 


trig 


2 


lam 


61/2 


0.245 


ct 


2 


lam 


77/2 


0.283 


sfcm 


2 


lam 


56 


0.383 


sm 


2 


lam 


58 


0.394 


fco 


2 


lam 


143/2 


0.484 


bee 


1 


com 


90 


0.567 


ct 


1 


com 


96 


0.594 


trie 


3 


lam 


337/3 


0.677 


trig 


1 


com 


115 


0.692 


sc 


1 


com 


128 


0.763 


hex 


1 


com 


138 


0.811 


bee 


1 


com 


169 


0.997 


ct 


1 


com 


180 


1.05 


sfcm 


1 


com 


210 


1.21 


fee 


1 


com 


229 


1.29 


ct 


1 


com 


243 


1.34 


hep 


2 


com 


263 


1.41 



